Relationship between non-exponentiality of relaxation and relaxation time at the glass 

transition 
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By analyzing the experimental data for various glass-forming liquids and polymers, we find that 
non-exponentiality (3 and the relaxation time r are uniquely related: log(r) is an approximately 
linear function of 1//3, followed by a crossover to a higher linear slope. We rationalize the observed 
relationship using a recently developed approach, in which the problem of the glass transition is 
discussed as the elasticity problem. 



Freezing a liquid to obtain glass may seem a familiar 
and conceptually simple process, and yet its theoretical 
description remains elusive. Such a description, as widely 
perceived, should provide a consistent theory for the two 
main properties that a liquid acquires in the glass trans- 
formation range: non-exponential relaxation and super- 
Arrhenius temperature dependence of relaxation time l| . 

At high temperature, a liquid under external pertur- 
bation relaxes exponentially fast: a relaxing quantity 
q{t) decays as exp(— (t/r)), where r is associated with 
the transition over a single activation barrier. This is 
known as Debye relaxation. On lowering the tempera- 
ture, relaxation changes markedly, and is described by 
a stretched-exponential function, q{t) oc exp(— (i/r)^), 
where < /3 < 1 P, 0| • The transition from Debye relax- 
ation to stretched-exponential relaxation (SER) marks 
the onset of glass transformation range. The trans- 
formation is complete when, by convention, the relax- 
ation time T increases to the experimental time scale 
of 100-1000 seconds, corresponding to the glass transi- 
tion temperature Tg. In the glass transformation range, 
T often increases faster than Arrhenius, and is well ap- 
proximated by the Vogel-Fulchcr-Tamman (VFT) law, 
T = Toexp(^/(r-To)) HI]. 

P and T are therefore two fundamental parameters that 
describe a liquid in the glass transformation range. A 
challenge for a theory of the glass transition is to pro- 
pose a description of these parameters. f3 and r have 
been discussed in a number of popular theoretical ap- 
proaches [l|, S S, 0, Hj @|- outstanding feature of 
these approaches is that P and r are often treated sepa- 
rately. One group of theories has offered the mechanism 
of the increase of r, and includes the Adam-Gibbs en- 
tropy theory 0], free volume theory [4], elastic models 
and other approaches (for a recent review, see Ref. [l|). 
Another group of theories has derived /3 for SER (see, 
e.g., Refs. 

In view of this, it remains unclear what the relationship 
between /3 and r is, or if one exists at all. On the other 
hand, because the increase of non-exponentiality and re- 
laxation time are the two signatures of glass transforma- 
tion, it is natural to ask if there exists a fundamental 



process of slowing down of molecular motion that affects 
both quantities. If a single mechanism affects both /3 and 
T, it should be reflected in a well-defined relationship be- 
tween these parameters at different temperatures in the 
glass transformation range. 

In this paper, we show that a universal relationship 
between (3 and r exists for all temperatures in the glass 
transformation range: log(T) is an approximately linear 
function of 1//3, followed by a crossover to a higher linear 
slope. We discuss the observed behaviour in the elastic 
picture of the glass transition. 

We have analyzed the experimental data on dielectric 
relaxation, includin g ou r recent results, as well as earlier 

data 0, Si, M MM, M,MM,M,M,MMM^ 



22|. At each temperature, (3 and r were determined from 



the location and width of the dielectric loss peak, respec- 
tively. Prompted by our recent work on glass transition 
Hi 24 1, we have plotted log(r) as a function of 1//3. In 



Figure 1, this dependence is shown for 15 different glass- 
forming systems. 

Figure 1 immediately highlights our central point, 
namely that a universal relationship /(/3,t) = can be 
identified for the studied systems in the entire range of 
glass transformation. First, at high temperature, log(T) 
is approximately proportional to 1/(3. Second, as the 
temperature is reduced, a crossover to another higher 
slope takes place. 

We note here that generally, presenting the data as in 
Figure 1 is attractive since it does not require scaling 
by Tg, an arbitrary quantity from a theoretical point of 
view. In addition to providing a universal relationship 
between /3 and r, the plot in Figure 1 can serve as a con- 
venient comparative representation of non-exponentiality 
in different systems approaching their glass transition. 

We also note that the universality of the pattern in 
Figure 1 is extended to the domain of high pressure. In 
the dielectric spectroscopy experiment, the increase of r 
at higher pressure can be counter-balanced by an increase 
of temperature. For many systems, it has been shown 25 1 
that different combinations of pressure and temperature 
that keep r constant, always give the same value of /3. 
This signifies the universal relationship between /3 and r. 
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FIG. 1: log(r) vs 1//9 for 15 glass-forming systems. Upper 
)anel: molecular liquids (■ o-terphenyl Q, • salol A PDE 
^, -A' 4-tert butyl pyridine 01 54 % chlorinated biphenyl 
lOl ]. , ♦ propylene carbonate dibutylphthalate [ll|], □ 

KDE 01, o BMMPC0), A 62% chlorinate biphenyl [H, V di- 
isobutylphthalate [Q). Lower panel: ▲ polyvinylacetate LZl, 
• polycyclohexylmethacrylate [l4| . I polyvinylethylene [15| . 
▼ bromopentane [l^ . Lines are guides for the eye. These 
represent the materials for which data are available covering 
a broad range and exhibiting "normal" behavior (i.e., /3 de- 
creases with decreasing temperature) . There are some glass- 
formers (e.g., NMEC [ilLMTHF 1,4-polyisoprene [H, 
polymethylphensiloxane |20l ]. 1-propanol ethanol 0) for 
which 13 is almost constant over the supercooled regime, and 
correspondingly there is a weak or absent crossover. In other 
materials, a neighboringdispersion, such as the normal mode 
(e.g., polyoxybutylene [2l|]) or a secondary peak (BMPC 
obfuscates accurate analysis of the a-process. For propylene 
glycol and glycerol [iJl, log(T) changes with (3, but does not 
show the linear behaviour as in this Figure. 



There is no a-priori reason why the location of the di- 
electric loss peak (r) should be correlated with its width 
(P). Hence the existence of a universal relationship be- 
tween (3 and r strongly suggests that they reflect the same 
slowing-down mechanism operative in the glass transfor- 
mation range. 

Figure 1 presents a challenge for a theory of the glass 
transition. In order to rationalize the observed be- 
haviour, one requires an approach to the glass transition 
that offers descriptions of both f3 and r. Some have been 
proposed 27|, and it would be interesting to see the 
predictions of these approaches regarding the observed 
behaviour. Here, we discuss the observed relationship 



based on the recently developed approach, in which the 
problem of the glass transition is discussed as the prob- 
lem of elasticity [2^, Below we briefly review the 
model, followed by a discussion of its prediction about 
the relationship between f3 and r. 

A glass is different from a liquid by virtue of its ability 
to support shear stresses on experimental time scales [l| . 
In our picture, we discuss that when considering stress re- 
laxation, a liquid can be treated as an elastic medium. In 
other words, we show that the problem of glass transition 
can be formulated as the elasticity problem. Important 
to this discussion is the elastic feed-forward interaction 
between local relaxation events (LREs) [23I 24 1. Lets 



consider the dynamics of LREs, induced in a liquid by an 
external perturbation (e.g. shear stress). Some time ago, 
Orowan's introduced terms of "concordant" and "discor- 
dant" events [2^: a concordant local rearrangement is 
accompanied by a strain agreeing in direction with the 
applied external stress, and reduces the energy and local 
stress. A discordant rearrangement, on the other hand, 
increases the energy and local stress. This has led to a 
general result that stress relaxation by earlier concordant 
events leads to the increase of stress on later relaxing re- 
gions in a system. Goldstein applied the same argument 
to a viscous liquid f29'| : consider a system under external 
stress which is counterbalanced by stresses supported by 
local regions. When a local rearrangement to a potential 
minimum, biased by the external stress, occurs (a con- 
cordant event) , this local region supports less stress after 
the event than before; therefore, other local regions in 
the system should support more stress after that event 
than before 

Let Ap be the increase of shear stress on a current 
LRE due to previous concordant LREs. If n is the cur- 
rent number of LREs, Ap is a monotonically increasing 
function of n. The increase of stress, on a currently re- 
laxing region increases its activation barrier V. It has 
been argued that V is given by the elastic shear energy 
of a surrounding liquid [sO]- This result was conflrmed 
by the experimental measurements of the shear modu- 
lus, showing that the activation barrier increases with the 
shear energy [sOl- Because, as discussed by Orowan and 
Goldstein, previous LREs reduce stress in the direction 
"concordant" to the external stress, the increase of shear 
stress on later rearranging regions consistently increases 
shear strain on them in the same direction, increasing 
shear energy and therefore V. The increase of V due to 
the additional stress Ap, AV, is given by work J Apdq. 
If Qa is the characteristic volume |3(ii, AV^ = Apqa, and 
we find V = Vo+QaAp, where Vq is the high-temperature 
activation barrier. 

Because Ap is a monotonically increasing function of n 
and V — Vo+QaAp, we find that V is also a monotonically 
increasing function of n. This provides the feed-forward 
interaction mechanism between LREs, in that activation 
barriers increase for later events. 
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It is important to discuss how the feed-forward inter- 
action mechanism operates on lowering the temperature. 
The elastic wave that propagates stress created by a LRE 
is stopped when another LRE takes place at the front of 
the wave. If r is the structural relaxation time, the elas- 
tic wave propagates without dissipation distance d — ct, 
where c is the speed of sound, d can therefore be called 
the liquid elasticity length. Because c is on the order 
of o/to, where a is the interatomic separation of about 
1 A and To the oscillation period, or inverse of Debye 
frequency (tq = 0.1 ps), 



d = a — 
To 



(1) 



Let drn be the distance between neighbouring LREs. 
dm is the distance between the centres of neighbour- 
ing molecular cages of about 10 A. At high temper- 
ature, when T « To, d < dm (see Eq. (1)). This 
means that neighbouring LREs do not elastically inter- 
act. Because events are independent, we obtain the ex- 
pected high-temperature result that relaxation is expo- 
nential in time and Arrhenius in temperature. Because 
T = To exp(y/fcT), a certain temperature always gives the 
opposite condition, d > dm, at which point the elastic 
feed-forward interaction mechanism between LREs be- 
comes operative. Note that the maximal time between 
two neighbouring LREs is given by t, hence d = ct < dm 
{d — CT > dm) also means that the time between the 
neighbouring LREs is shorter (longer) than the time of 
elastic propagation between the events dm/c. This is 
another way of showing that local events relax as inde- 
pendent at high temperature, but start to interact as the 
temperature is lowered. We have recently shown that this 
results in the crossover from exponential to stretched- 
exponential relaxation [23]. The crossover temperature 
Tc is calculated by putting d — dm in Eq. (1). 

To calculate how V depends on the current number 
of LREs n, we introduce the dynamic variable n{t), the 
current number of relaxing events induced by an exter- 
nal perturbation (i.e., in addition to thermally- induced 
events). We consider relaxation at constant tempera- 
ture. According to Eq. (1), this sets the range of the 
feed- forward interaction d. n{t) starts from zero and in- 
creases to its final value n^, n{t) ^ rir as t — > oo. Lets 
consider the current LRE to relax in the centre of the 
sphere of diameter d. As discussed above, all previous 
remote concordant LREs within distance d from the cen- 
tre participate in the feed-forward interaction, increasing 
stress Ap on the central region and hence increasing V for 
the central LRE. Ap can be calculated by integrating the 
contributions of remote concordant LREs. A straightfor- 
ward integration, together with V{n) = Vo + qaAp from 
the above, gives [23|: 




FIG. 2: /3 as a function of a, obtained from the solution of 
Eq. (3). 



where Vi — i: /2prqaApQd^ ln(2d/(io), do is on the order of 
the size of a relaxing region, is the density of relaxing 
regions, pr = Qur/Trd^, and Apo is the decrease of stress 
due to a remote concordant LRE. 

In Eq. (2), Vi depends on temperature through d (see 



Eq. (1)). As we have recently shown [^J, usmg t 
To exp(y/fcr) in Eq. (1) and eliminating d from Eq. (2), 
gives the VET law for t. Here, t corresponds to n = n^, 
i.e. maximal relaxation time of the system [2^. 

Eq. (2) describes the feed-forward interaction mech- 
anism in a liquid at T < T^. We are now set to write 
the equation that relates (3 and t. The rate of LREs, 
dn/dt, is proportional to the number of unrelaxed events, 
(n,. — n), and the event probability, p — exp(V/kT). Since 
V depends on n (see Eq. (2)), p becomes dependent on n. 
Introducing q — n/rir, and reduced time I/tq, we write: 



dt ^ 



q) cxp 



Vb + Viq 
kT 



(3) 



(2) 



Eq. (3) has two parameters, Vb/fcT and a = Vi/kT. 
We have recently shown 23[ that its solution is well 
approximated by the two-parameter SER, q(t) = 1 — 
exp (— (^/t)'^). Whereas t depends on both Vo/kT and 
a, P depends on a only; the smaller a, the larger (a = 
gives exponential relaxation, /3 — 1). We solve Eq. (3) for 
different values of a, fit the solution to the form of SER 
above, and find that = l + Ca, where C is a constant 
(see Figure 2). Joining this result with Vi oc ln(d) oc ln(T) 
(see Eqs. 1-2), we find that ln(T) is a quasi- linear func- 
tion of 1/(3. 

To discuss the crossover to the second higher slope in 
Figure 1, we note that our picture predicts crossovers of 
both T and (3 when at low temperature d reaches the 
system size L. The crossover of t takes place because 
V is proportional to VI oc ln(d) (see Eq. (2)) As 
long as d < L, V increases with temperature, because d 
is temperature-dependent (see Eg. (1)), giving the VET 
dependence of relaxation time [2J|. When d > L, V 
can not increase by way of increasing d, resulting in a 
weak temperature dependence of ^ at ci > L and thus 
the crossovers observed experimentally [1, 31 1 ■ The same 
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FIG. 3: Temperature dependence of in salol. Tl is the 
temperature of the crossover when d = L. 



reasoning is applied to the crossover of /3. Vi increases 
on lowering the temperature because Vi oc ln((i) in Eq. 
(2). This remains true as long as, on lowering the tem- 
perature, d < L. When d > L, Vi ln(L/(io), and is 
temperature-independent. Hence at d — L, Vi shows a 
kink and starts to saturate to a constant value. Because 
f3 decreases with Vi/kT (see Figure 2), we find that d = L 
should mark the crossover of f3 to the lower slope. Using 
salol as an example, d = L rs 1 mm in Eq. (1), and the 



VET parameters A = 839 K and Tq = 195 K for r [23, 
we find the crossover temperature of 247 K. This agrees 
well with the crossover of /? shown in Eigure 3. 

Because our picture predicts the crossover of /3 to the 
lower slope, the crossover of log(r) vs l/j3 to the higher 
slope also takes place, as is seen in Eigure 1. 

The crossover of /3 at cZ = L is expected to be more 
pronounced than that of r, because the activation bar- 
rier is the sum of a constant term and a temperature- 
dependent term hi^d/do) whereas f3 is solely defined 
hy a — Vi/kT (see Eqs. (2-3)). Consistent with this, we 
find that the crossover in Figure 1 is affected mostly by 
the crossover of (3. In Figure 3 we observe that for salol 
at low temperature, trails off at about 1.5, the same 
value of 1//3 associated with the crossover to the higher 
slope in Eigure 1. The crossover in Eigure 1 takes place 
in the temperature range in which r shows a crossover as 
well, although the latter is less pronounced as compared 
with the crossover of /3. 

In summary, we have shown that the non- 
exponentiality parameter /3 and relaxation time r are 
uniquely related in the glass transformation range, and 
rationalized this relationship in the approach that dis- 
cusses the problem glass transition from the elastic point 
of view. 

We are grateful to Prof. R. Bohmer and E. Roessler 
for providing their data in digital form, and to EPSRC 
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